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We perform a systematic analysis of several HBT parameters in heavy ion collisions from -Ebeam = 
2 AGeV to ^/snn = 200 GeV within the UrQMD transport approach and compare the results to 



ON 

experimental data where available. We find that the 'lifetime' of the emission source as calculated 

(N 



from r ~ y R"o ~ -^s 1 ' ^ s l ar g er than the experimentally observed values at all investigated energies. 

o 
o 



' The calculated volume of the pion source (Vj) is found to increase monotonously with increasing 



beam energy and the experimentally observed decrease of the measured Vf at AGS is not seen. 
Finally, we calculate the mean free path A/ = 0.5 — 1 fm of pions at freeze-out and find a good 
description of the experimental data above the AGS energy region, supporting the suggestion of a 
universal kinetic decoupling criterion up to the highest RHIC energies. 
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In order to create the theoretically predicted deconfined phase of Quantum Chromo- 
dynamics (QCD) heavy ions have been collided with energies from less than y/s ~ 2.5 GeV 
(SIS/FAIR energy regime), 2.5-20 GeV (FAIR/ AGS and SPS) up to 20-200 GeV (RHIC). 
Indeed, it seems that some nontrivial signals - such as charmonium suppression, relative 
strangeness enhancement, etc. - of the (phase) transition to the deconfined phase have 
been observed in heavy ion collisions (HICs) at SPS energies . Additional 

information about the matter created in such collisions can be obtained from the investiga- 



tion of the space-time structure of the particle emission source (the region o 



: homogeneity) . 



The established tool to extract this information is known as Femtoscopy 7|] or originally 



as Hanbury-Brown-Twiss interferometry (HBT) |3, Is , LlOf] - Experimentally this technique is 



quite often used to extract the information on the spatio-temporal evolution of the particle 



source, which has been scanned thoroughl y by se veral separate 



over the whole discussed energy region 
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351 ] . However, 



so far the measured excitation function of HBT parameters shows no obvious discontinuities 
within the large span of explored beam energies 0]. 

A comprehensive theoretical investigation on the excitation function of the HBT parame- 

Q 

ters is thus highly required but still absent so far [7[]. Recently, based on the Ultra- relativistic 
Quantum Molecular Dynamics (UrQMD, v2.2) transport model (employing hadronic and 
string degrees of freedom) (for details, the reader is referred to Refs. 23|, |24j, |25j, |26j) and the 



program CRAB (v3.0/3) 27|, |28|, |29J , we have investigated the transverse momentum, system- 
size, centrality, and rapidity dependence of the HBT parameters Rl, R q, R s (dubbed as 
HBT radii or Pratt radii), and the cross term Rql of pion source at AGS [30], SPS [3jJ and 



RHIC 



32| energies, respectively. In general, the calculations are satisfying and well in line 



with the experimental data although discrepancies are not negligible. Such as, I), the calcu- 
lated Rl and Rs values for Au+Au collisions at low AGS energies are visibly smaller than 
the data if the default UrQMD version 2.2 (cascade mode) is adopted. II), the HBT-'puzzle' 
with respect to the 'duration time' of the pion source, is present at all energies. In order 
to understand the origin of this HBT-'puzzle', some efforts were made, however, a complete 
understanding of this effect is still lacking. 

In this paper, we present pion interferometry results on the "duration-time" related quan- 



tity y Rq — R s , the freeze-out volume Vf, and derive the mean free path A/ of pions at 
freeze-out. The analysis is based on the comprehensive comparison of the excitation func- 



2 



tion of calculated HBT radii Rl, Ro, and Rs at small transverse momenta with data, in the 
FAIR/AGS, SPS and RHIC energy regime. The standard UrQMD v2.2 in cascade mode is 



employed throughout this paper to serve as a benchmark for further discussions 47|] . 



To calculate the two-particle correlator, the CRAB program is based on the formula: 

£/ k a = J d 4 x 1 d 4 x 2 g(x u pi)g(x 2 , P2) |0(q, r) | 2 
/ d i xig(xi, pi)/ d 4 x 2 g(x 2 , p 2 ) 

Here g(x, p) is the probability for emitting a particle with momentum p from the space- 
time point x = (r, i). 0(q, r) is the relative two-particle wave function with r being their 
relative position, q = P2 — Pi arid k = (p l + p 2 )/2 are the relative momentum and 
the average momentum of the two particles. Due to the underlying quantum statistics, 
this correlator is larger than unit at small q for bosons and can be fitted approximately 
by a Gaussian form. Using Pratt's three-dimensional convention (the LCMS system), the 
standard parametrisation of the correlation function in Gaussian form reads 

C(qo, qs, Ql) = 1 + Xex V (-R 2 L q 2 L - R 2 q 2 - R 2 s q 2 s - 2R 2 OL q q L ). (2) 

Here and Ri are the components of the pair momentum difference q and the homogeneity 
length (HBT radii) in the i direction, respectively. The pre-factor A is the incoherence 
parameter and lies between (complete coherence) and 1 (complete incoherence) in realistic 
HICs. The term R 2 OL is called cross-term and vanishes at mid-rapidity for symmetric systems, 



while it deviates from zero at large rapidities 31 
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34] 



We compare our calculations of the Pratt parameters of the pion source with experimental 
data for the following central collisions of heavy nuclei: 

1. Au+Au at the AGS beam energies E b = 2, 4, 6, and 8A GeV (< 11% of the total 
cross section cr T ), a rapidity cut \Y cm \ < 0.5 (Y cm = |log( ^ cm ^|| ), E cm and p\\ are the 
energy and longitudinal momentum of the pion meson in the center-of-mass system) 



is employed. The experimental (E895) data are taken from 12 ]. 



2. Au+Au at the AGS beam energy 11. 6A GeV (the < 5% most central collisions), a 
rapidity cut \Y cm \ < 0.5 is employed. The experimental (E802) data are taken from 



13|. 



3. Pb+Pb at the SPS beam energies E b = 20, 30, 40, 80, and 160A GeV (< 7.2%a T of 
most central collisions), apion-pair rapidity cut \Y W7r \ < 0.5 (Y^ = |log( f^^z^j^ ) 



is the pair rapidity with pion energies E\ and E2 and longitudinal momenta pp and 
P||2 in the center of mass system) is employed. The experimental (NA49) data are 
taken from 



4. Pb+Au at the SPS beam energies E b = 40, 80, and 160A GeV (the < 5% most central 
collisions), the pion-pair rapidity cut Y n7T = —0.25 ~ 0.25, —0.5 ~ 0, and —1.0 ~ —0.5 
are chosen. The experimental (CERES) data are taken from {if] ]. 

5. Au+Au at the RHIC nucleon-nucleon center-of-mass energies ^/snn = 30 (< 15%ar), 
62.4 (< 15%(7 T ), 130 (< 10%(7t), and 200 GeV (< 5%a T ). Here a pseudo-rapidity 
cut \r] cm \ < 0.5 (r] cm = ^\og(^^-) , (p is the momentum of the pion) is employed. The 



2 °\p— p 

experimental (PHOBOS, STAR, and PHENIX) data are taken from 
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2l|. 



Fig. [T] shows the excitation function of the calculated HBT radii Rl [in (a)], Rq [(b)], 



Rs [(c)], and the duration-time related quantity J Rq — R 2 S [(d)] at k T = 100 MeV (full 
lines, black) and 200 MeV (dotted lines, red). The experimental data within this transverse 

nn 

momentum region are shown for comparison. Since the experimental data from NA49 [14J,[15| 



and from CERES 16j collaborations overlap at beam energies 40, 80, and 160A GeV, we 
show the calculations and data with respect to CERES energies separately as dashed-dotted 
lines and open symbols. 

The calculated Rl (Fig. (U^a)) increases faster than R Q and Rs with increasing energies, 
which is also observed in data. Meanwhile, with increasing beam energies, the splitting of Rl 
with different becomes stronger, which can be attributed to the flow-dominated freeze-out 
scenario. This observation is also in line with previous results [3^], where it was found that 
the decrease of Rl with increasing at RHIC energies is stronger for central reactions than 
for peripheral ones. It should also be noted that the small kx behaviour of the correlation 
function is also affected by some other factors, such as the decay of resonances, potential 
interactions and/or the treatment of resonance life-times and widths in the medium. 

Fig. [1] (b) shows the calculations on Rq in comparison to the experimental data. Here we 
find that model calculation in data agree fairly well at AGS and SPS-NA49 but deviations 
from data are observable at RHIC energies. The comparison to the CERES data (with the 
appropriate cuts) however shows a rather large discrepancy between calculation and data. 



This deviations was also report in a previous study 3JJ and might hint to systematic differ- 



ences (apart from different experimental centralities and cuts) between the NA49 analysis 
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FIG. 1: (Color Online) Excitation function of HBT radii Rl [in (a)], Rq [(b)], Rs [(c)], and 



the quantity <J Rq — R? s [(d)]. The calculations are shown at k? = 100 ± 50 MeV (full line) 
and 200 ± 50 MeV/c (dotted line), respectively. The gray areas between the kx = 100 MeV 
and kx = 200 MeV lines are shown for better visibility. The data are at k? ~ 150MeV/c for 



n 



n 



reactions at E b = 2,4,6,8,4 GeV (AGS-E895) [Ij and 2J1 30, 40, 80, 160,4 GeV (SPS-NA49) [15j], 



at kx ~ 170MeV/c for reaction at. 



at E b = 11.6A GeV (AGS-E802) 



62.4, 200 GeV 



13, 



18 



20] 



SNN 



130 GeV 



2J, at k T ~ 200MeV/c for reactions 



13(], E h = 40,80,160,4 GeV (SPS-CERES) 



161 ]. and ^Jsnn 
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and the results measured by the CERES collaboration. The source of this difference remains 
therefore unclear. Recently, the NA57 collaboration published the Kt dependence of the 
HBT radii (with Kt up to 1.2 GeV/c) in Pb+Pb collisions at 40A GeV, and it was found 



35] 



that the NA57-data, especially the Ro, are in line with the NA49-data 

In Fig. [11(c), we notice that the calculated sideward radii R$ are in qualitative agreement 
with the data but seem to be 15% smaller for almost all energies. Furthermore, the increase 



of the measured Rs at low AGS energies can not be reproduced [30] and might be due to 
the omission of potential interactions that gain importance at low beam energies. Detailed 
investigations will be presented in a future publication. The excitation function of the HBT 
radii from several systems (from light to heavy) inspired by the NA49-future collaboration 
4j might also provide new insights into this problem, and predictions are in progress. 



The HBT duration time "puzzle", i.e. the fact of the theoretical quantity J R Q — R$ 
being larger than extracted from the data, is present at all investigated energies (see Fig. 



CD (d)): The calculated values of J Rq ~ R% are about 3.5 ~ 5 fm while the measured ones 



are 



30 



.5 ~ 4 



31 



m. Many efforts have been put forward over the 



36, 



37 



38 



ast years to clarify this issue 



31] , we suggested that a 



39j, E.g., in our previous works [30|, 
sizeable amount of interactions between particles at the early stage of the reaction either on 
the mean field and/or (non-perturbative) partonic level seem to be relevant to understand 
the phenomenon. In fact, all straightforward cascade transport approaches (also those with 
partonic interactions with perturbative QCD cross sections) fail to describe the quantity 
Rq — Rg, as well as the elliptic flow [41], over the whole energy range. 
Fig. [2] (a) shows the excitation function of the pion source volume Vf at freeze-out, 



calculated as 



42 1 



V f = (2n)lR L R 2 s . (3) 

Note that the radius Ro is not considered to calculate the pion freeze-out volume since it 
contains the contribution of the temporal extent of the pion source. Fig. [2] (a) shows clearly 
that the UrQMD cascade calculations do provide a reasonable freeze-out volume for the pion 
source at RHIC energies. At SPS energies, the agreement is fine with CERES data while it 
slightly underpredicts those of NA49. Towards even lower energies, the model underpredicts 
the measured freeze-out volume due to the omission of the strong interaction potential and 
other in-medium effects. E.g. at E b = 2A GeV, the measured Vf is about 2 — 3 times larger 



than calculated value. As studied in 



301 ]. a mass-dependent lifetime of resonances accounts 




FIG. 2: (Color Online) (a): Excitation function of the pion freeze-out volume Vf (according to Eq. 
[3]) at transverse momenta between = 100 MeV and 200 MeV (gray area), compared with data in 
this /^-region, (b): Excitation function of the mean free path Xf of pions at freeze-out (according 
to Eq. (31) at the same transverse momenta. The experimental value for A/ at ^snn = 200 GeV 



is obtained with the help of recent dN/dy data in 



from 



431 ] . at all other energies the Xf data are taken 



for an improvement of the HBT-radii at small kx and hence reproduce the data better. 
From the discussions above, it is clear that the major part of this difference is related to the 
sideward radius Rs- 

We are now ready to estimate the mean free path Xf of the pions at freeze-out from the 



following expression 



42 1 



Vi 



f 



No Nnv-xn + N^a-n 



(4) 
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with the averaged pion-nucleon cross section <jntt = 72 mb and the averaged pion-pion cross 
section a n7T = 13 mb (Note that within the present model calculations these values are 
slightly energy dependent. However, here we have adopted the explicit numbers from Ref. 



42j to compare to the results presented there). The nucleon and pion multiplicities and 



AL are calculated as 



dy 
and 

a r /7T~ dNpi ons 

N n = y th ■ v 27r 



N N = y th ■ V2tt — \ VnM , (5) 



(6) 

Uni id 



dy 

using the assumption of a thermal equilibrated system at freeze-out with a temperature 
Tf = 120 MeV. Here, yth is the estimated thermal homogeneity scale in rapidity at a certain 
far and Tf, and is given by the expressions: y t u = arctanh^/?^)), with {(3 t h) = \Jl + {l) 2 / (j) 
and (7) = 1 + 1/3 (Ki(m T /Tf)/ K 2 (rriT/Tf) — 1)+Tj /m T . Here K n {z) is the modified Bessel 



function of order n and my = ^va\ + k\. At = 100 MeV and 200 MeV, the calculated 
homogeneity lengths in rapidity y t h are 0.98 and 0.81, respectively. dN / dy\ Vmid is the rapidity 
density of pion (nucleons) at mid-rapidity [48[. Recent calculations using the present UrQMD 
transport model [25[ , have shown that the calculated pion and nucleon yields are reasonably 
in agreement with data. 

Fig. [2] (b) shows the excitation function of Xf of pions at freeze-out. It is seen that the 
theoretical Xf value increases gradually from ~ 0.5 to ~ 1 fm from AGS to highest RHIC 
energies with a weak dependence on kx- The experimental values of A/ are also between 
0.5 — 1 fm. The calculated Xf below the AGS energy deviates from the data due to the 
smaller calculated sideward radii (resulting in a too small freeze-out volume). 

Certainly, the assumption of a constant kinetic freeze-out temperature Tf = 120 MeV in 
the calculations of y t t over the whole energy range is not justified, especially for the collisions 
at low beam energies. By inserting an energy dependent temperature (assuming a kinetic 



freeze-out temperature of 70% of the chemical freeze-out temperature shown in 44j) into 
Eqs. |5] and El one finds that at lower AGS energies, the number of nucleons and pions in 
the pion source volume is visibly reduced due to the decreased Tf. As a result Xf increases 
and a flatter distribution of Xf as a function of the beam energy is obtained in the present 
calculation. 
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The observation (both experimentally and theoretically) of a nearly energy independent 
mean free path on the order of 0.7 fm at pion freeze-out is rather surprising. Physically it 



has been interpreted as a rather large opaqueness of the pion source at break-up 45|, |46 ] . 
To summarize, we survey the excitation function of the HBT radii Rl, Rq-, and R$, 



the quantity J Rq — R s , the volume V/, and the mean free path A/ of pions at freeze-out 
for heavy systems with energies ranging from lowest AGS to the highest RHIC energies. 
Generally, the model calculations with UrQMD v2.2 (cascade mode) are in line with the 
data over the whole inspected energy range. Although discrepancies especially in the lower 
AGS energy region are found and have to be resolved. Especially the HBT parameter R$ and 
the HBT duration-time related "puzzle" deserve further attention. Finally we re- interpret 
the measured and calculated radii in terms of a mean free path for pions and kinetic freeze- 
out. Here we find a nearly constant mean free path for pions on the order of Xf = 0.7 fm 
indicating a significant opaqueness of the source. 
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We have found [30] the treatment of the mass-dependent resonance lifetimes can better de- 
scribe the HBT radii for AGS energies (at small transverse momenta), meanwhile, the HBT 
time-related puzzle can be better understood with the consideration of a potential interaction. 
These effects are, however, not taken into account in this paper for the sake of clarity. 
[48] In Ref. [l2], the experimental N nuc i eons is set to two times of total amount of emitted proton 
and anti-proton yields and N p i ons is three times of negatively charged pion number. In our 
calculations, N nuc i eons represents the sum of protons and neutrons, as well as anti-nucleons. 
Note that this introduces a small difference into the analysis due to iso-spin effects, especially 
at very low energies. 
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